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Abstract 

We show that a bipartite state on a tensor product of two matrix algebras 
is almost surely entangled if its rank is not greater than that of one of its 
reduced density matrices. 

1 Introduction 

1.1 Background 

Recently, Arveson ]2| considered the question of when a bipartite mixed state of 
rank r is almost surely entangled, and showed that this holds when r < d/2 where 
d is the dimension of the smaller space. In this note we show that this result holds 
if r < d, with d now the dimension of the larger space. 

* Partially supported by National Science Foundation under Grant DMS-0604900. 
'''Partially supported by an NSF grant, a FRG-NSF grant and a BSF grant 
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We will use results from [11J on entanglement breaking channels and exploit 
the well-known isomorphism between bipartite states and completely positive (CP) 
mapsQ We will first consider states associated with completely positive trace- 
preserving (CPT maps) and then find that extension to arbitrary bipartite states is 
quite straightforward. 

If the rank of a bipartite state ^ab is strictly smaller than that of either of its 
reduced density matrices, then the state must be entangled. This is an immediate 
consequence of well-known results on entanglement, and seems to have first appeared 
explicitly in [12]. We include a proof in Appendix |A] for completeness. This allows 
us to restrict attention to the case in which the ranks of the reduced density matrices 
are equal, with one of full rank. 

Although it seems natural to expect that this result is optimal, recent results of 
Walgate and Scott [19] suggest otherwise. Let the Hilbert spaces TCa and TCb have 
dimensions (1a and d B respectively. It follows from a result proved independently by 
Wallach [20] and by Parthasarathy [15] for multi-partite entanglement that when 
s > (d,A — 1)(g?b — 1) any subspace of Hab = 7~La ® 7~Lb with dimension s contains 
some product states, and that this bound is best possible, i.e., if s < {(1a — 1)(g?b — 1) 
then there is some subspace of dimension s with no product states. 

Walgate and Scott extended this by proving [19], Corollary 3.5] that if a subspace 
of Ha <S> 7~Cb has dimension s < (gU — 1)(cIb — 1) then, almost surely, it contains 
no product states. For a bipartite state jab with rank r < (c?a — 1)(^b ~ 1), it 
follows that range of jab almost surely contains no product states, which implies 
that a bipartite state jab with rank r < (d,A — 1)(c?b ~ 1) is almost surely entangled. 
Alternatively, one could apply [19], Theorem 3.4] directly to ker(7^#) to reach the 
same conclusion. 

When cLa > d B > 2, this result is stronger than ours, but for a pair of qubits, 
d-A = ds = 2 our result is stronger. Moreover, it is easy to extend our qubit results 
to the general case of bipartite states with rank r = dA > ds > 2, providing a proof 
quite different from that in [15] . Although our measure is constructed differently 
from that used in [2], our approach is similar in the sense that we show that in a 
natural parameterization of the set of density matrices, the separable ones lie in a 
space of smaller dimension. 

In the next half of this section, we review relevant terminology, and describe the 
notation and conventions we will use. Qubit channels and states are considered in 
Section |2j and the general case in Section [3] We conclude with some remarks about 
other approaches, and the question of the largest rank for which the separable states 

1 This isomorphism is usually attributed to Jamiolkowski [13 or to Choi [7], who used it to 
characterize the complete positive maps on finite dimensional algebras. However, it seems to 
have been known to operator algebraists earlier and appeared implicitly in Arveson's proof of 
Lemma 1.2.6 in [T] 
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have measure zero. 



1.2 Basics and notation 

In this paper, we consider maps $ : B(Ha) B(TCb) and identify them with 
bipartite states or, equivalently, density matrices in B(Ha) <S> B(Hb) via the Choi- 
Jamiolkowski isomorphism as described below. Our primary interest is the situation 
in which Ha = Cd A , m which case we can identify B{Cd) with M d , the space of 
d x d matrices. However, we will also have occasion to consider either Hilbert space 
H as a proper subspace of for some d. 

We will identify a state with a density matrix, i.e., a positive semi-definite op- 
erator p with Trp = 1, in B(H). To an operator algebraist this corresponds to 
the positive linear functional on the algebra B(H) which takes A i— > Trp A. In the 
physics and quantum information literature, a density matrix (or, more properly, a 
density operator) is often referred to as a (mixed) state on H (because the density 
operator acts on H. ) 

When Ha = Cd A and Hb = Cd B , we write <3> : M dA i— > M dB . In this case, 
let {ej} and {f m } denote orthonormal bases for C dA and Cd B respectively. The 
isomorphism between states and matrices arises from the fact that 

Tr\f m )(f n \m ej ){e k \) (1) 

can be interpreted as either 

(i) the matrix representative of the linear map $ : M^ A i— > Md B in the bases 
\fm)(fn\ and |ej)(e fe | for M dB and M djv respectively, or, 

(ii) the density matrix jab of a state on Cd A <8> C^ s with elements [7Afi]jm,fcn i n 
the product basis \ej <S> fm)- 

Conversely, any state on Cd A <8> Cd s defines a CP map. We describe this well- 
known fact in detail in order to establish some conventions for interpretations of 'Ja 
and 7 B . Observe that (ii) is equivalent to writing ^ab as a block matrix of the form 

Iab = ^Y^\ej)(e k \®P jh = ^^|e,)(e fe | ® $(| ej )(e fe |) (2) 

jk jk 

with the block Pj k = $(|ej)(efc|) the matrix in Md B given by the image <&(|e_j)(efc|). 
One can write an arbitrary matrix in Md A <g> Md B in the block form Yljk \ e j)( e k\ Pjk 
and then define <3>(|ej)(efc|) = Pjk and extend by linearity or, equivalently, 

^{A) = Y,^kP jk (3) 

jk 
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when A = Y,jk a jk\ e j)( e k\- 
Observe that 

1b = iTr A7AB = ^^$(|e fc )(e fe |) = ^$(/ A ) (4a) 

k 

1A = ^Tr B7AB = ^^|e J )(e fc |Tr$(|e J )(e fc |) (4b) 

jk 

and that this implies the following: 

a) $ is unital, i.e., $(I A ) = Ibi if and only if 1b = o~Ib, and 

b) $ is trace-preserving (TP), i.e., Tr B $(X) = Tr A X V X e 13 (Ha), if and 
only if 7 A = ±I A . 

When Mrf or 13(H) is equipped with the Hilbert-Schmidt inner product, one can 
define the adjoint, or dual, of a map $. We denote this by $ and observe that this 
is equivalent to 

TrB^(A) = Tr [$(5)] f A (5) 

A matrix $ is TP if and only if its adjoint $ is unital. 

It is a consequence of Theorem 5 in [7] that the extreme point^] of the convex 
set of CP maps for which 7 A = $(/b) = p have a state representative (often called 
the Choi matrix) with rank < rank p. We prefer to consider CPT maps and regard 
the density matrices with rank < d A as an extension of the set of extreme points. 
As shown in Appendix |Bj this corresponds to the closure of the set of of extreme 
points. We let T>c denote the set of density matrices in B(Hc) or M dc and T>c(r) 
to denote the subset of rank r. We also define the following subsets of T>ab(t)- 

V A (p; r, s) = {7ab e V A b ■ rank7 AB = r, rank7 A = s and ^a = p} (6a) 
Va(t, s) = {'Jab e £>ab : rank 7^ = r and rank7 A = s} (6b) 

Although the sets in ^ above are subsets of V A b C B(Ha) ®B(Hb) — M dA ®M dB 
we use the subscript A to emphasize that we impose conditions only on the marginal 
7 A . When rank pi = rank p 2 = d A , the map 

Iab i-> (pl /2 Pi 1/2 ® Ib)iab (Pi l,2 p l 2 2 ® Ib) (7) 

gives an isomorphism from V A (pi] r , d A ) to Va(P2, r, d A ) and each of these is isomor- 
phic to Va(j^Ia', dA, dA) which is isomorphic to the set of CPT maps $ whose Choi 
matrix has rank d A - We will let «S A (p;r, s), etc. denote the corresponding subsets 
of separable state in ([6]). 

It will be useful to introduce the notation for the map that takes a density 
matrix p 1— > T^pT. 

2 Choi's condition for true extreme points is implicit in Theorem 1.4.6 of [T]. 
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2 Maps with qubit inputs 



2.1 Canonical form and parameterization 



Now consider the case of CPT maps on qubits for which TCa = TCb = C 2 . As 
observed in [2], these maps can be written using the Bloch sphere representation 
in the form 



$ (w I + 22 w kVk) = w I + y^(t fc wp + \ k w k )a k 



(8) 



where cr& denote the three Pauli matrices. Necessary and sufficient conditions on 
tk, Xk which ensure that $ is CP are given in [IS] . The form ^ is equivalent 
to representing $ by a matrix T with elements tjk = |Tr (Tj$(crfc) so that, with 
subscripts j, k = 0, 1, 2, 3 and the convention J 2 = do 

/l 0\ 

ti Ai 

t 2 A 2 

\t 2 A 3 / 



(9) 



As shown in [TH Appendix B] an arbitrary unital map on qubits can be reduced 
to this form by applying a variant of the singular value decomposition to the 3x3 
submatrix with j,k G {1,2,3} using only real orthogonal rotations. Given the 
isomorphism between rotations and 2x2 unitary matrices, this corresponds to 
making a change of basis on the input and output spaces TCa = Cd A = C 2 and 
TCb = Cd B = C 2 respectively. Thus, for an arbitrary unital CP map $ one can 
find unitary U, V such that T^t o $ o has the form (J8|) or, equivalently, a matrix 
representative of the form (|9|. 

It was shown in [18] that the maps with Choi rank < 2 are precisely those for 
which the form ^ becomes 



/ 



T 



1 






ysin u sin v 





cos ii 

cos v 











(10) 



cos u cos v J 



withj^j-u, v in (— it, tt] x [0, it]. Moreover, as shown in [IB], the entanglement breaking 
(EB) maps are precisely the channels which have either cosw = or cosv = 0. 



It follows from (10) that every element of V A {\I] % 2) can be represented by 
a triple ((w, v), U, V^jconsisting of a point in R 2 , and two unitary matrices U, V. 



3 The interval for u is shifted from that in 18J. However, the interval [0, 7r] for v was incorrectly 
stated as [0,7r) in [15] , 



5 



However, some care must be taken so that each element of 7^(1/; 2, 2) is counted 
exactly once. It suffices to restrict (u, v) to the rectangle 



-r- r 7T-, r 7T-, 

A = [0, o] x [0, -] 



'11} 



2 J 2 ■ 

Suitable rotations will give all allowed negative values of the non-zero elements in 



(10), as well as even permutations of and Problems with overcounting occur 



only on the lines u = 0,v = 0,u = v. To deal with this we define 
A = {(u,v) : < u <^,0 < v <^,u ^ v}. 



(12) 



(The line segments on the boundary with u — f and v — f are included in A as 



shown in Figure 2.1 



(o,s: 



(0,0) 





( 7T TV \ 




" ✓ 

1 ✓ 
1 ✓ 

1 ✓ 

1 ✓ 

. ✓ 

1 ✓ 

1 ✓ 

1 ✓ 

* ✓ 

1 ✓ 

1 ✓ 

, ✓ 

1 ✓ 

1 ✓ 

1 ✓ 

1 y 




k. 





• u 



Figure 1: The rectangle A corresponds to the shaded region. The dashed lines are 
not in A. The lines u = | and v — | correspond to the EB channels. 

Because different pairs of matrices U, V may give the same channel on the 



lines not included in (12), we define equivalence classes as follows. Let 7Z t (with 
t = x,y, z) denote the subset of SU{2) corresponding to the rotations around the 
indicated axis. We write (U, V) ~ (U' , V) if there is an R t G TZ t such that U' = R t U 
and V = R t V or, equivalently U'W = V'V^ G TZ t , and denote the quotient space 
(SU(2) x SU(2))/7Zt. With this notation, we now make some observations 

a) The subset of EB channels consist of those channels for which either u = | or 
v = -■ 

b) The line u = v corresponds to the amplitude damping channels (It is well- 
known that only the case u — v = | is EB; this is a completely noisy channel 
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mapping to a fixed pure state.) From (10) one sees that these channels are 
invariant under rotations about the z-axis, and the set of amplitude damping 
channels in V A (\l]2,2) is isomorphic to (u,u) x (SU(2) x SU{2))/K Z . 

c) The line segments with u = and v = correspond to phase-damping chan- 



nels. From (10) one sees that these channels are invariant under rotations 
about the x and y-ax.es respectively. Thus, the set of phase damping channels 
in Va(\I') 2, 2) is isomorphic to 

{(0,v) : v e (0, f ,]} x (SU(2) x SU{2))/U X 

U {(u,0) :ue(0, §]} x (SU(2) x SU(2))/K y . 

d) The point u = i> = gives the identity channel, for which rank 'Jab = 1- 
Thus Va(^I', 2, 2) is isomorphic to 

A x STf(2) x 517(2) |J {(«, w)} ue(0)f ] x (ST/(2) x SU(2))/K Z 

|J {(u, 0)} ue(0lf ] x (ST/(2) x SU(2))/K y (13) 
|J{(M}*=(o,§] x (^(2) x SU{2))/K X 
and, 5a(|/; 2, 2), the subset of EB channels in Va{\I', 2, 2), is isomorphic to 
{(u, f ) : u e (0, |)} x SU(2) x £17(2) 

U {(I > v ) :w e (°. I )} x 5f/ ( 2 ) x 5f/ ( 2 ) U (°> \) x ( 5f/ ( 2 ) x SU ( 2 ))/ n * 
|J (|, 0) x (ST/(2) x ST/(2))/7^ |J (|, |) x (SU(2) x SU(2))/H M . (14) 

2.2 Construction of a measure 

Let m2 be the normalized Lebegue measure on A and v 2 the normalized Haar 
measure on SU(2). Then the product measure Jl = 7712 x i/ 2 x i/ 2 defines a probability 
measure on Q 2 — A x SU(2) x SU(2). Although every point in Q 2 corresponds to 
an element in Va{^}\ 2,2), it can happen, as described above, that more than one 
point corresponds to the same CPT map $. Therefore, to define a measure on 
Va (|7; 2, 2) we use the map g : fi 2 — > Va(^I] 2, 2) which takes 

((u,v),U,v) h-Tyo^oTrVt (15) 



where denotes the CPT map whose Choi matrix is given by (10). The map g 
is surjective which allows us to define a measure /x on all sets X C V A (\l]2,2) for 
which g~ x (X) is measurable by 

fi(X) = Jl(g- l (X)). (16) 
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Since g is surjective, g 1 (P^(|J; 2, 2)) = Q 2 which implies that Va{\I\ 2, 2) is mea- 
surable and ii{Va{\I] 2, 2) = 1. Thus, /i is a probability measure on Va(\I\ 2, 2) 

Moreover, the entanglement breaking channels satisfy 
M(5 A (|/;2,2)) 

= //({(u, |) : u G [0, |]} x 517(2) x 517(2) U {(§ ,u) : u G [0, f } x 5(7(2) x 5*7(2)) 
= 0-l-l + 0-l-l = (17) 

Thus, we have proved the following 

Theorem 1 A CPT map $ : M 2 i— ► M2 0/ Choi-rank 2 is almost surely not EB, 
or, equivalently, a state jab on C 2 <8> C2 which has rank 2 and 7 a = \l is almost 
surely entangled. 

Since, the unitary conjugations have Choi matrices of rank 1, and correspond 
to the set (0, 0) x 5*7(2) which has measure zero, we have also proved the following 
result, which we state for completeness. 

Theorem 2 A CPT map $ : M 2 hh> M 2 of Choi-rank < 2 is almost surely not EB, 
or, equivalently, a state jab on C 2 <8> C 2 which has rank < 2 and ja = \l is almost 
surely entangled. 



2.3 Removing the TP restriction 

We would like to extend the results of the previous section to 

Theorem 3 If a state 'Jab on C 2 (g> C 2 has rank 2 and ja a lso has rank 2, then 
'Jab is almost surely entangled. 

Proof: As observed after (JTl), V a^PY^Aa) — T- > A(P2', r ,dA); Indeed, the CP maps 
corresponding to states in T>a(p; 2, 2) have the form <£>oT v / 2 ^ with $ CPT, although it 

might seem more natural to consider the dual T ' r^p o $ which takes 1 1— > dAP- Next, 
observe that any density matrix p G M^ A of rank 2, can be written as U ( g 1 ^_ x ) U' 
with x G (0, |) and *7 G SU(2); the case x = | gives |7 independent of U. Thus 
the set of density matrices p G M^ A of rank 2 is isomorphic^] to 

\l U (0,|) x 5*7(2) (18) 



^Hcre we use the fact that a x pa x exchanges the eigenvalues. This is quite different from the 



situation in ( 12 1 where we could not assume u < v because the permutation in S3 which exchanges 



1 <-> 2 can not be implemented with a rotation. 
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and the set of bipartite density matrices Va{2, 2) (for which rank jab = rank 7^ = 2) 
is isomorphic to 

P a (±I;2,2)|JPa(|/;2,2) x (0, \) x SU{2). (19) 

To define a measure on this set, let mi denote normalized Lebesque measure on 
(0, \) and let X 2jt be defined using product measure so that 

A m= f^ x ™i x ^)W XeV A &;2,2)x(0,l)xSU(2) 
2A ' [(l-t)fi(X) XeV A (ll;2,2) 1 ' 



where we can pick any t G (0, 1] and /i is the measure defined in Section 2.2. Then 
the subset of EB channels Sa{2, 2) has measure 

\ 2>t (S A (2,2)) = fi(SA(ll;2,2))+fi(SA(ll;2,2)) mi (0,l)u 2 (SU(2)) 

= t • + - • 1 • 1 = QED (21) 

independent of t G (0,1]). We can drop the requirement that ja has rank 2 by 
observing that extension to all jab of rank 2 requires only that one replaces (0, |) 
on the right side of (19) by [0, |). Thus, we can conclude that 

Corollary 4 If a state jab on C 2 <S> C2 has rank 2, then jab is almost surely 
entangled. 



2.4 Two-dimensional subspaces of C^. 

We can use the isomorphism between C 2 and any Hilbert space of dimension 2 
to replace either Ha or by a two dimensional subspace of C^. However, for 
later use, we now want to extend our qubit results to the somewhat more general 
situation of the set of all CPT maps $ : C2 1— > Cd B whose range has the form 
£>(span { |v 1), |f2)}) with \v\), \v 2 ) G Cd B - Here, we do not fix the range, but consider 
all CPT maps whose range corresponds to some two-dimensional subspace of Cd B - 

Observe that in the polar decomposition Tyt o $ o Tjj leading to the canonical 



form rt8j) we need only replace V by an isometry V : C 2 1— » He- Then in (13) and 



(14), the first use of SU(2) in each subset must be replaced by Vd which is defined as 
the subset of d x 2 matrices satisfying V'V = I 2 . By Theorem A. 2 of [2], Vd can be 
given the structure of a real analytic manifold with a probaility measure Vd (which 
is unique if it is required to be left-invariant under SU(d)). Although V is not a 
group, we can define equivalence classes as before with(V, U) ~ (V, U') if there 
is a R t G !Z t such that V = VR t and U' = UR t . Then the previous arguments 



go through with SU(2) x SU(2) replaced by V d x SU(2) in Section 2.1 and the 



corresponding use of v 2 in Section 2.2 by Vd- 
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3 General maps 



3.1 CPT maps with d A > 2. 



We now assume d& > d% > 2 and extend these results to bipartite states on 



J d A 



Cd B with 7^4 = j^Ia- We begin by considering a CPT map $ : M t 



M, 



with Choi-rank dA- By Theorem 5C of [TT], which is equivalent to Corollary 14, $ 
can always be written in the form 



(22) 



where [gk] is an orthonormal basis for Cd A , but the states ipk G Cd B need not be 
orthogonal or even linearly independent. In the basis gk, the Choi matrix for $ has 
the form 



lAB= dkYl \9k)(9k\ ® |^fc>(^fc| 



(23) 



which implies that 7^ is block diagonal with each block a dn x <^b rank one pro- 
jection. Let us first assume that ipi and ip2 are linearly independent. 

Now let Pk = IV'&X'W and write (23) explicitly in block form, as 



1AB 



(Pi 








. 


.. \ 





P2 





. 


.. 








p. 


. 


.. 


v° 







. 


• • pd A ) 



(24) 



and consider a density matrix of the form 

(Q 




dA 





P3 
















\ 



p d A J 



(25) 



where Q G M 2 ®M^ fl is a positive semi-definite 2<ie x 2<ie matrix of rank 2 satisfying 
I2. Now a density matrix of the form (25) is separable if and only if |Q is 



Tr B Q 

separable. However, \Q is a density matrix of the form considered in Section 



2.4 



Let yd A ({gk}, {4>k}) denote the subset of Va( 



1 I a, dA-, d/ 



dA 



consisting of density 



matrices of the form (25) or, equivalently, 

dA 

y dA ({9k},{^k}) = {Q © Yl \9k)(9k\ ® iV'fcXV'fcl = Q e ^b.'&bQ = ^2} (26) 



fc=3 
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where © denotes the direct sum and 

X AB = B(spaa{|0i>, \g 2 )} ® £(span{|</>i>, \ip 2 }}). (27) 

The set of projections {ip^i^kl G M dg is isomorphic to S 2 d B -i, the £ 2 unit sphere 
in H 2dB . For a given \gi), \g 2 ), the set depends only on span {(V'l)) IV^}} and 
not the choice of individual vectors. Therefore, we can identify each point in 

Q dA = A 2 x V dfl x SU(2) x SU{d A )/SU{2) x x . . . x 

V v ' 

= A 2 x V dB x 5E/(dA) x 5 Mfl _i x . .. x S 2dB -i (28) 

S v ' 

d A -1 

with a density matrix ^ AB in = [j {gk} ^ k] ^d A ({Pfc}, {?M), the set of all density 
matrices of the form (25). (Note that S 2 d B -i occurs d A — 2 times in (28) correspond- 
ing to the choices of if)k for k = 3,4 . . . n. The set ({(?*:}, {ipk}) depends only on 
span {|^i), |^2)} — range V with V G V dB , with non-orthogonal vectors |?/>i), |^ 2 ) 
associated with non-unital qubit channels via isomorphism.) 



Let m 2 and v d be measures as in Sections 2J2 and 2.4 let ^ be normalized Haar 
measure on SU(d) and let n 2dg _ 1 be a probability measure on S 2 d B -i- We define a 
normalized measure /x on by the product measure 

Ji = m 2 x v dB x z/ 2 x z/ dA/2 x n Ms _i x . . . x n 2dB -i (29) 



d A - 2 



To obtain a measure on y dA we proceed as in Section 2.2. Let G : Vt dA 1— > 3^2 A be 
the map that sends an element ((u, v), V, U, ■ ■ ■ , \i/>d A )) to the corresponding 
density matrix in y dA and define 

MX)= A 7(G- 1 (X)) (30) 

whenever X C y dA for which G -1 (X) is measurable. 

As explained above, Corollary 14 implies that S A (-^I A ;d A ,d A ) C y dA . Then, 



proceeding as in ( |17j ), one finds 

»(S A (±I A ; d A , d A )) = • 1 • 1 • l d ^ 2 = 0. (31) 

Moreover, for any reasonable extension of /x from y dA to all of V A (-^-I A , ; d^, d^), 
the EB subset will still have measure zero. In particular, one could simply let 

«X) = h X) ' f ' Y C ^ , „ (32) 

\o if x cv A (j- A h]d A) d A )\y dA 

and note that u is absolutely continuous with respect to any other extension of /x. 
Thus, we have reduced the general case to that of d A = 2 and conclude that 
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Theorem 5 Let ^ AB be a state on C dA ® C dB which has rank oIa > d B > 2 and for 
which 7a = j^Ia- Then jab is almost surely entangled. 

Remark: The assumption that ip\ and ip2 are linearly independent can be 
dropped because that case corresponds to u = v = | in ([8j and is included implicitly 
in our analysis. The set of channels for which all are identical also has measure 
zero, except for the excluded situation d B = 1, for which all states are separable. 



3.2 Reduction of the general case to CPT 

As observed earlier, when rank p = d A 

V A (p;dA,d A ) = {{^d A p® I B )lAB{\/dAP® Ib) ■ Jab e V A (-^I A ;d A ,d A )} (33) 

is isomorphic to V A (j^I A ;d A ,d A ). But parameterizing the set of density matrices 
of rank d A is a bit more subtle than for d A = 2 because of the need to consider 
degenerate eigenvalues, for situations beyond j-I. However, this only affects a set 
of measure zero and can be dealt with as in the preceding sections. To describe the 
set of density matrix of rank d A consider the set of vectors 

Z = {z = (Ci,C 2 , ■ ■ -Q A ) : < Ci < C 2 < • • • < Q A ,^2Ck = l} (34) 

k 

in the positive facet of the i\ unit ball of Rd A - We can associate each z G Z with 
a diagonal matrix A z so that that the map h : (z, U) i— > UKJJ^ takes Z x SU(d A ) 
onto V A (d A ), the set of density matrices in M dA with full rank d A . Since we can 
identify Z with a subset of H dA -i, we put normalized Lebesque measure m dA _ 1 on 
Z, and let 

VdA (X) = (m dA . 1 xu dA )(h-\X)) (35) 



whenever X C X'a(gU) and /i _1 (X) is measurable. Then it follows from (31) that 
for any extension u of p, the product measure to x r\ dA gives a measure on Va^Ai <Ia) 
for which the separable states Sa^a^a) have measure 0-1 = 0. Thus, we have 
proved 

Theorem 6 If a state jab on C dA ®C dB has rank dA > d B > 2 and rank(jA) = dA 
then jab is almost surely entangled. 



3.3 Further results 



Theorem 11 states that if the rank of ja is dA and the rank of jab is strictly smaller 
than dA, then jab is entangled. Thus r < dA implies that V A (p] f, d A ) consists 
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entirely of entangled states. If we combine this with our results for r = s = we 
obtain several additional theorems, which we state for completeness. 



Theorem 7 Assume > ds > 2. If a state jab on M dA <8> M dg has rank jab < 
dA = rank ja, then 'Jab is almost surely entangled. 



By using the isomorphism between Cd and any Hilbert space of dimension d we 
can restate this by letting Ha = range ja and Hb = range 7# and considering jab 
as a state on B(Ha) <8> 13 (Ha)- 

Theorem 8 If a state ^ab on Md A (8> Md B has rank ^ab < rank 7^4 and rank 7^4 > 
rank 7^ > 2, then ^ab is almost surely entangled. 

We also find that we can eliminate the need to consider the rank of 74. 

Theorem 9 Assume dA > dg > 2. If a state ^ab °n Mj A ® Md B /ias ran/c 7ab < 
dA, then ^ab is almost surely entangled. 



Proof: Let Z denote the closure of (34). Since this simply replaces the strict 
inequality < Ci by < 0, the set Z x SU(d) includes all density matrices in Md A 
so that 



Z\Z x SU(2) = h- 1 ({p e V A ■ rank p < d A }) 



(36) 



Now extend the measure rj in (35) to all of V>a- The set of all separable states jab 
with rank jab < dA is Sa^a) = U s <d A <Sd A (dA, s). The subset of separable states 
with rank 7^4 < satisfies 

[J <Sd A (d A , s) C {p E V A ■ rank p < d A } (37) 

s<d A 



But 



Thus 



Vd A {{p e V A : rank p < d A }) = va dA -x(Z\Z) v dA (SU(2)) = 0-1. (3* 



{u dA x 7] dA )(S(d A )) = (uj 4a x Vd A )(<S dA (d A ,d A )) + x r) dA )( (J S^cZ^s 

s<c2a 

QED (39) 



< 0-1 + 1-0 = 
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4 Final comments 



4.1 Remarks on measure 

If we apply the argument used in to prove Theorem [9] to the subset of states with 
7^4 = j^Id A or equivalently, combine Theorems |5j and jllj we obtain the following 
result which we state in terms of channels. 

Corollary 10 Assume oIa > ds > 2. Then the set of CPT maps $ : Cd A i— > Cd B 

whose Choi matrix has rank r < oIa is almost surely entanglement breaking. 

As shown in Appendix B, the closure of the set of extreme points of CPT maps 
$ : Cd A i— > Cd B is precisely the set of channels whose Choi matrix has rank < oIa- 
Because the extreme points of a convex set lie on the boundary, their closure always 
has measure zero. Thus, Corollay 10 is a special case of a well known, more general 
fact from convex geometry. An alternative, and somewhat simpler, approach to 



proving Theorem 5, would be to use this fact together with Theorem 15 However, 
we feel that it is useful to see the specific paramaterizations which lead to our 
results. In our approach, one sees that everything really follows from the basic 
paramaterization of extreme points for qubit channels, and the fact that (up to sets 
of measure zero) the relevant sets of bipartite states can be parameterized as direct 
products on which we can put product measures. 



One could extend Corollary 10 to the set of CP maps for which $(-/#) = p 
with p G Va{t) fixed, again using the fact that the closure of the set of extreme 
points has measure zero. Then we can conclude that the subset of separable states 
U s < r S(p; r, s) has measure zero with respect to a measure on U s < r V(p; r, s). How- 
ever, we can not go directly from this observation to Theorem [9] by taking the 
Up6©A because the set T>a is uncountable. One would still need the argument in 



Section |3.2| What this observation about extreme points does tell us is that our 
results are not sensitive to the choice of measure. The fundamental issue is that the 
bipartite states can be parameterized as a smooth manifold on which the separable 
ones correspond to a space of smaller dimensions. 

There is one unsatisfying aspect of using the the inverse image to define a mea- 



sure, as in (16); namely, that it does not reflect the fact that different unitaries give 



the same map on some lines in A. An alternative would be to first define sepa 

7T 

' 2 



rate measures on the different regions in (13), e.g., on {(«,«)} , n tt, x (SU(2) x 

SU(2))/1Z Z , use the product measure m\ x v z where mi is normalized Lebesque 
measure on (0, |y and v z is Haar measure on the group (SU(2) x SU{2)) /1Z Z . One 



could then combine the measures on the four subsets in (13) as in (20) using, say, 
weights 1 — t x — t y — t z , t z , t x , t y with t m > and Y2m=i — !• However, given that 
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each of the line segments with u — f , u — v, and v — | have measure zero in A, 
the most natural choice weight would be t m = 0, equivalent to simply omitting the 
corresponding channels (or states). 



In fact, all regions which a quotient space is needed, as in Sections 2.2 and 3.2 
have measure zero in our inverse image approach. Intuitively, one would like to 
simply observe that we can identify T>a(\I] 2, 2) with a subset of [0, |] x [0, |] that 
satisfies A C T>a(^I', 2,2) C A and then observe that since 

//(A) <//(D4§/;2,2)) <//(A). (40) 

and yu(A) = /i(A) = 1, one must have /i(P^(|/; 2, 2) = 1. But to use this approach, 
one must establish that Pa(|-^;2,2) can be identified with a measurable subset of 
A. 



4.2 Optimality 

It is natural to ask if the results in Theorems [7] and [8] are optimal. For > 2, 
it is clear that the results which follow from those of Walgate and Scott [19J are 
better. Thus, the question becomes whether or not rank jab < (^a — 1)(^b — 1) 
is optimal. This does not follow from the subspace theorems in [19] because when 
rank ^ab = 2 > ((1a — 1) {ds — 1) the product states can form a set of measure zero in 
a subspace of Ha ® W^. However, we know that the separable ball in B(Ha ® 7~Cb) 
has strictly positive measure [U EJ [21] so that the optimal rank must be strictly 
smaller than dAds- 

In the case of qubits, we know that Theorem[3]is stronger than the results implied 
by Walgate and Scott [19J, and that when rank jab = 4, the separable states have 
strictly positive measure. If we restrict attention to those states jab with rank 3 
and 7a = 7b = \l or, equivalently, the unital CPT maps with Choi-rank 3, we can 
use the familiar picture of a tetrahedron [HI [161 HH] • The rank 3 states correspond 
to the faces, and the subset of separable states on each face to the smaller triangle 



whose vertices are midpoints of the edges as shown in Figure 4.2 Thus, the unital 
CPT maps with Choi-rank 3 have measure 0.25 with respect to all the unital CPT 
maps on qubits. However, it is open whether or not this holds when the restriction 
to unital maps is removed. Thus, the question of whether Va(^I'i 3, 2) has measure 
zero or positive measure seems to be open. 

Acknowledgment: It is a pleasure to thank Jonathon Walgate for bringing 
[19J to our attention, to Michael Wolf for providing reference p2] and to Michael 
Nathanson for assistance with figures. 
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Figure 2: The left figure shows the tetrahedron of unital qubit channels with the 
octahedron the EB subset. The right figure shows one of the faces of the tetrahedron, 
corresponding to channels with Choi-rank 3, with the shaded region the subset of 
EB channels 

A Some separability theorems 

For completeness, we now state and sketch proofs of some results that are well 
known and/or proved in [TTj . The first result appeared as [T2l Theorem 1] in a 
slightly stronger form. 

Theorem 11 If rank jab < d,A — rank ja, then jab is not separable. 
Proof: First observe that jab is separable if and only if 

Iab = j~ a {Ia ® 1aY 1/2 1ab(Ia ® 1aY 1/2 (41) 

is separable. But ja = ^~Ia- Now both the reduction and majorization criteria 
[6j [9] for separability of a state pab imply that the largest eigenvalue must satisfy 
||Pab||oo < 1 1 PB I |oo- But rank jab = rank jab < oIa implies that 7ab has at least 
one eigenvalue > Thus ||7ab||oo > j- — 1 1 1b 1 1 00, an d it follows that both jab 
and jab are entangled. QED 

When rank 7^ < oIa, one can regard the underlying Hilbert space as TCa to be 
range 7^4 = (ker r yA) ± - One then obtains 

Corollary 12 If rank jab < rank 74, then jab is not separable. 

The following Lemma goes back at least to \lU\ and a simpler proof was given 
in [11]. To emphasize that one need not assume oIa = &b (and because of typos in 
[TT] ) we include a full proof here. 
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Lemma 13 Let pab be a density matrix on TIa <8> Kg. If pab is separable, pab 
has rank d, and pa has rank d, then pab can be written as a convex combination of 
products of pure states using at most d products. 



Proof: Since pab is separable it can be written in the form 

k 



Pab = y^Aj \ai)(a,i\ g> \bi)(bi\. 



(42) 



i=i 



with ||aj|| = = 1. Assume that k > d and that pab can not be written in 
the form (42) using less than k products. Since pA has exactly rank d, there is 



no loss of generality in assuming that the vectors above have been chosen so that 
|ai), \a 2 ), ■ ■ ■ \a d ) are linearly independent. Moreover, since pab has rank d < k, the 
first d+1 vectors |aj) <S> \bi) must be linearly dependent so that one can find ctj such 
that 



d+1 

Now let {|efc)} be an orthonormal basis for Tis- Then 



(43) 



d+1 



(44) 



Since the first d vectors |o 3 -) are linearly independent, the solution of Y^j x j\ a j) = is 
unique up to a multiplicative constant. Applying this to the coefficients in (44) one 
finds that there are numbers i/& such that aj(ek, \bj) = v^Xj. Let \u) = ^2 k Vk\ek)- 
Xj\u). Since multiplying Xj by c, changes Uk — > -Vk, one can assume 
v\\ = 1 = Then <Xj\bj) = Xje^lu), and 

v). Therefore, , one can rewrite (42) as 



Then aj\bj) 

that Xj has been chosen so that 



ctj 7^ implies \bj) 



Pab= X i \ a M a i\ ® \ b j)( b i\ + X i\ a M a 3\®\ u )( u \ 



(45) 



j:aj=0 



Suppose that t of the otj are non-zero. Since the vectors {|aj) : a.j ^ 0} are linearly 
dependent, the density matrix ^2j. a j L \j \aj)(aj\ has rank strictly < t and can be 
rewritten in the form J2l=i X 'j\ a 'j) ( a 'j\ usm g onr y s <t vectors \a'j) . Substituting 



this in (45) gives pab as linear combination of products using strictly less than k 



contradicting the assumption that (42) used the minimum number. QED 
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Corollary 14 If ^ab is separable and 7^4 = j^Ia, then ^ab can be written in the 
form 

iab = ito*)<0*l ® I^X^fcl ( 46 ) 

k 

with gk an orthonormal basis for Cd A 

Proof: Since jab is separable it is a convex combination of projections onto product 
states and can be written in the form 

Iab = ^ £fclfffc ® tpk) (9k ® V'fcl ■ (47) 

Since rank 7^ is (i^ by assumption, it follows from Lemma [A] that we can assume 
k — 1,2 . . .dA (duplicating terms if < dA are needed). But then, the assumption 

j I I A = lA = Y,^\9k)(gk\ (48) 

k 

holds if and only if £ fc = i V A; and the vectors g% are orthonormal. QED 



B Closure of the set of extreme points 

It is often useful to consider the set of all CPT maps with Choi rank < dA- In [18J 
these were called "generalized extreme points" and shown to be equivalent to the 
closure of the set of extreme points for qubit maps. This is true in general]^] We 
repeat here an argument form [T7]. Let £(o?a, ds) denote the extreme points of the 
convex set of CPT maps from Md A to Md B . 



Theorem 15 The closure S(dA, ds) of the set of extreme points of CPT maps $ : 
Md A 1— > Md B is precisely the set of such maps with Choi rank at most dA- 

Proof: Let A\~ be the Choi-Kraus operators for a map $ : M^ A 1— > Md B with Choi 
rank r < dA which is not extreme, and let be the Choi-Kraus operators for a 
true extreme point with Choi-rank dA- When r < d& extend A\~ by letting A m = 
for m = r + l,r + 2, . . . dA and define Cfc(e) = A^ + tB^. There is a number e* such 
that the d\ matrices Cj(e)Cfc(e) are linear independent for < e < e*. To see this, 
for each Cj(e)Cfc(e) "stack" the columns to give a vector of length d\ and let M(e) 
denote the d\ x d\ matrix formed with these vectors as columns. Then det M(e) is 



5 Arveson [3] has pointed out that Theorem 15 can also be proved using results in [2J. 
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a polynomial of degree d\, which has at most d\ distinct roots. Since the matrcies 
AjA k were assumed to be linearly dependent, one of these roots is 0; it suffices to 
take e* the next largest root (or +1 if no roots are positive). Thus, the operators 
Cj(e)Cfc(e) are linearly independent for e G (0, e*). The map p i— > J2k Ck{t)pC k {e) 



3 

is CP, with 



C\{e) C k (e) = (1 + e 2 )I + e(A[B k + B\A k ) = S(e). 

k 

For sufficiently small e the operator S(e) is positive semi-definite and invertible, 
and the map $ e (p) = C fc (e)5(e)- 1 / 2 pS(e)- l ' 2 Cl(e) is a CPT map with Kraus op- 
erators C k {e)S{e)^ 1 / 2 . Thus, one can find e c such that e G (0, e c ) implies that 
$ e G £(d A ,d B )- It then follows from lfm $ e = $ that $ G £(d A ,d B )- QED 
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